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BUEGSaS APPROHtMATXON FOR TWO-DIMENSIONAL 
FLOW PAST AN ELLIPSE 

By 

J. Mark Dorrepaal* 

INTRODUCTION 

ECforts during this- phase of the research have concentrated on analys- 
ing the Burgers flow: past a circular cylinder, coopering predictions of the 
Burgers model with those of Oseen flow past a circle, and deriving exact 
solutions of the-flow equations- for. elliptic geometries. The results are 
provided in the following sections of this report. The first section fol- 
lowing the List of Symbols describes a motivation-for studying Burgers flow. 
Then "Solutionl* outlines a solution technique which works equally well for 
Oseen or Burgers flow past a circular cylinder.- This ia followed-by sec 7 
tiohs which describe the separation. behind the cylinder, the drag experi- 
enced by the cylinder, and asymptotic behavior far from the cylinder. The 
section titled "Burgers Flow and Oseen Flow: How Do They Differ?" shows 

that the predictions of Burgers flow near the cylinder provide a substantial 

inq>rovement over those of Oseen flowu Finally, "Burgers Flow Past an El- 

* 

liptic Cylinder" gives the formulation and solution of the equations of 
motion f >r flow past an ellipse. 

LIST OF SYMBOLS 

a cylinder radius 

A(6) coefficient in least squares fit 

Cq drag coefficient 

E coefficient in vorticity expansion 

n 


*.issociate Professor, Department of Mathematical Sciences, Old Dominion 
University, Norfolk, Virginia 23508. 
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(C,n) 

(W,X) 


radial aigenfunctloh raauUing £tom aeparacioa of 
variables 

modified MBChieu function 

Green's function . 
unit vector in a*4itection 
modified Bessel function 
pressure 

rear stagnation, point of the cylinder 
polar coordinates 
Reynolds number 

Critical minimum Reynolds number 

Reynolds number basod on Che semimajor axis of the 
ellipse 

Mathieu function 

angular eigenfunction resulting from separation of 
variables 

free*>8tream velocity 
fluid velocity 
Cartesian coordinates 

Reynolds number exponent in least squares fit 
coefficient in infinite linear system of equations 
Kroneckar delta 

delta function 

gradient operator 
Laplacian operator 

elliptic coordinates for flow parallel to major axis 
elj^iptic coordinates for flow perpendicular to major 
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V kindioafcio visoooicy 

p fluid denaicy 


♦ (r»®) 
♦h<r) 
X(R) 
X(M,M 

'Ke.n) 

0 ) (r,0) 


atreaa cofliponenfia 

harmonic coAjugata o£ potontial flow aCream function 

coefficient of higher order term in aaynptotlc 
expansion of stream function 

coefficient of second term in asyfi^totic ei^ansion of 
stream function 

stream function for potential, flow past an ellipse 
perpendicular to its major axis 

stream function for Oseen or Burgers flow past a circle 

stream function for Burgers flow past an ellipse 

stream function for potential Jlow past an ellipse 
parallel to its major axis 

vorticity 


MOTIVATION 

The ttondimensional Mavier**Stokes equation has t he form 

R(v • ^)v ■ -^p + V^ V (1) 


where R, v» p are Reynolds number, fluid velocity, and pressure, 
respectively* The velocity vector nmist also satisfy the continuity 
equation 

^ • V ■ 0 (2) 

which is guaranteed if we use the following representation: 

V • curl {^(r,9)k} (3) 
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where t|»(r,0) le the stream functiont By substituting equation (3) Into 
equation (1) and ellmiitatlng p(r,0), we can write the Navler-Stokes equa 
tlon In the equivalent scalar £orm 


(I) m t|> 


(4) 



(5) 


dr 36 


„ 3i|> 3 

iTil 


This suggests an. iterative process 


72 li • (I) 

k 


j^V2 + - j W 


k+1 


- 0 


k • 0 j l» 2|^* ^ ♦ . 


( 6 ) 

(7) 


which defines a sequence {t|»|j(r,0)} of spatially uniform ap.proximations to 
the solution of equation (1). It is assumed that under appropriate Condi- 
tions this eequenc ill converge to the solution of equation (1). 

To begin the iterations we must propose initial values of wo(r,0) 

and t|»o(r,®) consistent with equation (6). Since we are considering 
uniform flow past a circular cylinder r ■ I, the flow is irrotational at 

infinity and it is reasonable to choose wo(r,6) s 0. If we choose t|»o “ 

r sin corresponding to a uniform stream^ the problem for (wi» 70 i® 

Oseen flow. On the other hand, if we choose 7o “ 

ponding to continuous inviscid flow past the cylinder, the resulting problem 
for (wj, ,|,i) is Burgers flow. Both flows should behave similarly at infi- 
nity. But in Oseen flow vortlcity l« convected through the cylinder while 
in Burgers flow it is convected around the cylinder. This suggests a 
difference in the two flows near the circle. 


4 


SOLUtXON 


With'k M fl. and (jig ■ t ain Qf aquafctott (7) bec owei the Oeeon vorticlty 
equ.^tieAi 

fv* 4 . RiilLlL." a 6 oi 0 L-lw(r, 0 ) " 0 (d) 

L r 30 3rJ 

The Burgers vertioity equation has the form 



j ain-0 L- - r| 


1 cos 0 1- u(r,6) ■ 0 

(9) 

L Vr rS, 

f 3d ' 

^ r^ / 

■ 3r J 



Both solutions are of the form 


QO 

o)(r.e) - I K F„(r) T„(0) 

n-1 « “ « 


( 10 ) 


where E are constants to be determined and (^, P^, T are given in 

n n n 

In table 1. 

Table 1. functions occurring in vorticity expression. 


Method 

<Kr,0) 

F„(r) 

T^C9) 


Oseen 

1 Rr ^os 0 


sin n 6 


Burgers 

R cosh 2 cos 0 

Cekj^ (a» - i a^) 

*®n -ia2) 

a ■ In 1 




The £unecions Qelt^ and ee^ are Mathleu £unationi| the fiormer behaving 
asymptotically like exp <-R cosh a) as a « in-r-de^a, the latter being odd 
and periodic in 6 with period 2t. 

The stream function corresponding to w(r|6) satisfies the Poisson 
equation (4) and tho boundary conditions 

t|»(l,a) -li (1,0) - 0 (ID 

dr 

A dreen's function defined by 

V* G(r,0}r^,0^) ■ - ^ 0 (r - r^) 0(0 - 0^) 

O(r,0;l,0 ) - 0 
o 

is used to solve for <|>(r,0) with the result: 

'i'<V®o^ • ** ®o C w(r,0) 

•if A 

• G(r,0jr ,0 )r dr d0 (14) 

O 0 


( 12 ) 

(13) 


Where 


G(r.0;r^,0 ) - - i. in . 

4it 

To obtain the coefficients 
condition in equation (14): 


1^2 4. y2 „ 2 rr„ cos (^ ~ 1 

o o ^1 

(rr + I - 2 rr cos (0 - 8 ) 

0 O o J 

in equation (10), we invoke the no-slip 


0 


(l,0_) ■ 2 sin 0 

0 


♦ f r «(r,0) 
-r ^ 


3G 

dr 


(r,0;l,0^) r dr d0 


(15) 
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Expanding ^ in a Fourier iine eerioi in 0^ and equating the resulting 
*'o 

coefificienta^ o£ ain k6 (k ■ 1,2,3,.,.) on the right side of equation (15) 

0 

to aero, wo obtain an infinite linear systeia of the form 


00 - 


I 

n"l 


E r, 

n kn. 


kl 


k * 1,2,3,... 


(16) 


where 


r . ^ f* r t‘"'‘ *' °°* ° 

7 0 ‘ 




sin n6 


• sin k0 dr d0 (Oseen) 


r -1 f” r* ^ Gek.. 


z r r 

^ ft a 0 \ 4 / 


• 8 


in k0 se^j ^0 * * — da.de (Burgers) 


If equation (16) is truncated at k • n 8, good results are obtained for 
Reynolds numbers in the range 0 < R < 4. 

SEPARAtlOM 


Both Oseen and Burgers flows exhibit separation on the downstream side 
of the cylinder provided the Reynolds number exceeds some critical minimum 
^c. Ibid value is found by solving the equation 


9(i> 

— (l.0;Rj - 0 
30 ® 


(17) 


Yamada (ref, 1) has shown that R^ ■ 1.51 for Oseen flow and our calcula- 
tions verify this result. For Burgers flow we find R^. • 1 . 12 . A numeri- 
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cal value obtained by Underwood (ref# 2) from the Navier''8Col(ea equation ia 
Rc * 2.8b. 

Xt ia expected that the Burgera reault would be Idea than the numerical 
value* The oonvective velocity field in Burgera flow ia continuoua inviacid 
potential flow peat the cylinder, and-thia violatoa the no-alip condition at 
the cylindoela aurfaco. The velocity field which aolvoa the full Naviec- 
Stokea equation aatiafioa thia condition. Thua, convection offeeta near the 
cylinder are more dominant in Burgora flow than in. Navior-Stokoa flow, and 
any phenomena related to convection, auch^ aa separation, should occur at 
lower Reynolds numbers in Burgers flow* 

The fact that the Burgera result is less than the Oseen value is a 
little surprising, but can be explained. Separation, begins at the rear 
stagnation-point P of— the cylinder where locally the flow appears as in. 
figure 1. At the onset of separation a bubble of circulating fluid forms 
about E. The direction of motion along the axis of symmetry inside the 
bubble is opposed to that outside the bubble (fig, 2). In Oseen flow the 
convective velocity field is constant in magnitude and perpendicular to the 
cylinder boundary in the vicinity of P as shown in figure 3, Oseen con- 
vection therefore will deter the establishment of reverse flow at P be- 
cause it directly opposes the direction of fluid motion along the axis of 
symmetry inside the separation bubble. In contrast, the convective velocity 
field in Burgers flow is parallel to the cylinder boundary near P; in 
fact, its magnitude vanishes at the point P Itself (fig. 4). Burgers 
convection does not oppose the establishment of a separation bubble about P 
in the way Oseen convection does, and we would therefore eiqpect separation 
to initiate in Burgers flow at a lower Reynolds number. 

As R increases beyond R^ ■ i.ij Burgers flow, the wake grows in 
siae. When R ■ 2.0 the flow appears as in figure 5. The length of the wake 
is PQ ■ 0.53 where OP is the unit of length and <80P ■ 34.8*. Point t, 
whose B -coordinate is 83*, marks the location where the fluid pressure along 
the boundary is a minimum. The flow from T to 8 is against an adverse 
pressure gradient. 
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A circular cylinder in a uniiorm etreaa experiencea a force in Che 
direction of the flow at- -infinity. The magnitude of this force ia obtained 
by ;kntograting the component Of Che atreac veotor^n^ the direction of the 
uniform atroam about tho circumforeneo ef-.tho circle. The drag ia thoroforo 


D 



coa 9 


0 - ain e) de 
« /r-l 


(18) 


where 0 ^ are the atreea components and p, v, U are fluid 

density, kinematio. viaoosityi and free-stream velocity, respect ivoly* The 
stress cofi^onents can be calculated from equation (14) with the result 


2pvU 0.k„(0. 


fir R cos 9 
• / e 

o 




sin 6 d6 


(W 


The drag coefficient is defined to be 


C 


U 


(p U2 a) 


( 20 ) 


where a is cylinder radius. Table 2 provides drag Coefficients as compute 
ed from equation (19) for Reynolds numbers in the range 1 < R < 2. 


Table 2. Drag coefficients for 1 < R < 2. 


R 

1.0 1.12 1.25 1.50 2.0 

Cjj 7.76 7.29 6.86 6.22 5.35 
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These vslyes agree well wich Che experiiriSncs o£ TrlCCon (see refis. 3 end 4). 
Table 2 will be expanded shorcly Co include Reynolds numbers in Che range 
I < R < 4 and drag opeg£ieienCs £or ellipcic eylinders. 

ASYMPTOTIC FOfttiOR-RLOW-RAR FROM CYLXNU6K 


The s«4uCion o£ Che full Navior-SCokos equation past.a .inito obstacle 
predicts o boundary layer sureounding the obstacle in which the flow is 
rotational and tho- velocity gradients are Large. Outside this boundary 
layer the flow is essentially irrotational. Mathematically this means Che 
stream function outside the boundary layer is harmonic. Since Oseen flow 
and Burgers flow- are spatially uniform linear models of a Navier-Stokes 
flowi we would expect them both to exhibit this behavior* 

If the Gr-eenLs function 0(r|8ir ,6 ) in equation (14) is etpatv^ne in 

o o 

a Fourier series t the coefficients of cos kS (k ■ 1,2,...) ''U voirUh^ ieav~ 
ing the stream function in the following form: 


r 

<|>(r^,0o) • “ — I. •” ® r^ ^ ai(r,9) sin kS dr d0 


ff k-l k ° ° o 


I I i r sin k6^ r ! ° w(r,d) 


V k*l k 
• sin k0 dr d0 


o I 


( 21 ) 


This expression can in turn be written 


♦ Bin 0^ - sin 0^ 

0 0 O 0 o 


1 jt f o 

ir 0 


r // w(r,0) sin 6 dr d9 ♦ I 


-I 


( 22 ) 

(cont'd) 
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-i I i-lnke. 
ir k*2 k 


jv I'^Q yl~k k8 dr d8 

0 ^ 






I 

k 


8ia-k8^ 


r 

Jit J 0 ^l+k gin Jj0 dr d6 

0 ^ 


( 22 ) 

(concl'd) 


Vhe. three expressions in brsekets in equation (22) all have finxte limits as 

'o * ** Taking this limit, therefore, will yield the second temuin the 
asyng>totic ex pan sion 6f the stream function far from the cylinder* The 
result is 


w 

r|»(r_^,0^) ~ r sin 9. - — x(R) I ~ sin k0 


o o o 


° It .k-1 k ° 


•"(y 


(23) 


Now it is known that 


OD 

I i sin k$ ■ (it - 8 ) 0 < 0. < It 


k*l k ° 2 ' ° 


(24) 


- (-It - 0 ) -It < 0^ < 0 
2 ° ® 


(25) 


Thus the stream function has the asymptotic form 


♦^'o»®o^ ~ ’'o • 


in 0^ - i x(R) 


(26) 


with the plus sign being taken if 0 < < it and the minus sign if -it < 

0 <0. The second term in equstion (26) can be made analytic along 0 ■ tt, 

o o 

but suffers a jump discontinuity across Oq ■ 0. 
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BKpAnsiona oil Che •Cream £an«tion fer 6o ” ^ Cq Hnice reveal 
no dlaconClnuiCyi howcvert The dlaconCinulCy in equaCion (26) ia a propnrCy 
o£ Che aaympeocio expanaion of i(i(roi ^uC not -oX ^he acream 

funcclon iCaelX. 

Thin rachor aCrange behavior can be explained. Prom equnCion (10) che 
vorCiciCy can be ahowi Co behave eaympcocically like 

-1/2 Rr (I - C08 8) 

w(r,8) f(6) aa r « (27) 

>F 


where €(0) e o, f'(0) ^ 0. Thus Che vorticiCy decays exponencially as 
r • provided 0 0. Along 0 ■ 0, however , we have u(r,0) ■ 0 for ail r 

and U (r,0) ■ 0 (r"^^*) as r Thus (r,0) decays algebraically along 

0 ■ 0. If r ■ li » I, Che funcCion (L,0) will be exponencially small as— - 


0 *** 0, buC will jump Co a much larger value when 0 • 0 since ics decay alo ng 
Chis ray is so much- slower-. This behavior- is noc unlike a delca function, 
and, since Che aCream funcCion is relaCed Co Che vorCiciCy Chrough Poisson's 
equaCion, ic is noC surprising ChaC a sCep function behavior appears in Che 
asympCoCic expansion of t|)(r,6) along 0 ■ 0.~ It must be stressed, however, 
ChaC Chis is a feature of Che asymptotics only. The stream funcCion itself 
is analytic in the fluid domain. 


The next term in equation (23) can be calculated and is of Che form 

♦j(R)r^"^ sin 0^. In fact, all subsequent terms are of the form ♦j^(R)r~*^ 
sin k0g (k ■■ 1,2,3,...). The asymptotic expansion of the stream function 

is therefore harmonic (except possibly along 0q ■ 0) and gives us the 
potential flow far from the cylinder. Table 3 contains some computed exam> 
pies. 


BURGERS PLOW AND OSBBN FLOW: HOW DO THEY DIPPER? 

The previous section shows that one distinct advantage of a Spatially 
uniform linear approximation to the Navier-Stokes solution is the abi;ity to 
calculate from it the resulting potential flow at Infinity. This potential 
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Tabid* 3* Aayi6pt@fclc expanalona Coc latga 


Oaeen Flow 
R ■ 1,51 



1 - 0.5995 

R • 2.0 

2.9247 


0.5777 

R - 3.0 

2.4888 


0.5485 

R • 4.0 

2.2481 


0.5302 

R - 5.0 

2.0918 


0.5171 


Burgers Flow 



flow cannot be obtained a priori without aome knowledge of the flow near the 
cylinder because the flow near the cylinder determines the outer edge of the 
boundary layer, ^ich in turn defines the region of irrotstional flow. 
Without an approximation like Oeeen flow or Burgers flow, therefore, the 
problem of finding <t>(r,6) far from the cylinder reduces to finding a harmo- 
nic function in a region whoso boundary is unknown. 

Roving obtained the asymptotic expansion given in the previous section, 
however, it is possible to define the locotion of the outer edge of the 
boundary layer as the curve along which this asymptotic expansion vanishes. 
This amounts to finding a boundary which the potential flow does not pene- 
trate. Slippage along this boundary is, of. course, permitted. The boundary 
so defined determines the displacement body which the potential flow far 
from the cylinder "seesj’ The-displacement body includes the cylinder, its 
wake, and the boundary layer surrounding the cylinder. 

By setting the asymptotic expansions given in table 3 equal to sero and 
solving foe r, we -can obtain, approximations to-ldte displacement bodies for 
the various flows. A typical example is -given in figure 6. Note that the 
displacement body is semi- infinite with its thickness at infinity being 
twice the value of the (il ** 4) ** coefficient. Since the cylinder boundary 
is given by r ■■ 1, the thickness of the boundary layer is easily calculated. 
Tables V-and 5 compile these results for a variety of Reynolds numbers and 
locations, along the cylinder boundary. Note that 6 « tt refers to the for- 
ward stagnation point on the cylinder boundary. 

The last Column in tables 4 and S is a least-squares fit of the data 
given in each row. In table 4 the value of a depends on 0. In table S, 
however, _the value of a hovers about ■ 1/2 regardless of 9. This 
suggests that the boundary layer thickness is inversely proportional to R^^^ 
in Burgers flow but not in Oseen flow. 

Studies of the nonlinear boundary layer Clow past a semi- infinite flat 
plate show unquestionably that the boundary layer thickness behaves like 
Our work obtains a similar result for a bluff body using a linear 
model in which vorticity iS convected around the obstacle's boundary. The 
agreement between this pred’*tion of Burgers flow and that of nonlinear 
analysis regarding boundary layer behavior indicates that the Burgers model 
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0.356a 0.2360 


0.36A9 0.2427 

0.3904 0.2637 

0.4363 0.30L4 

0.5086 0«3610 

0.6185 0.4515 

0.7858 0.5895 

1.0483 0.8066 

1.4840 1.1688 

2.2783 1.8334 

3.9903 3.2763 

9.4140 7.8764 


0.1691 

0.1254 

0.1750 

0.1308 

0.1936 

0.1478 

O.2270 

0.1784 

0.2796 

0.2267 

0.3597 

0.3001 

0.4818 

0.4121 

0.6743 

0.5887 

0.9962 

0.8846 

1.5892 

1.4314 

2.8834 

2.6288 

7.0282 

6.4777 


0.7984/R1*^336 

0.803d/R^*“*’ 

0.8230/R^‘®®‘*® 

0.8647/R®*^^*^ 

0.9415/R®*®’®® 

1.0706/R®*’®®® 

1 . 2796 / 8 ®*^®**** 

1.6182/R®*®®®* 

2.1876/R®*®®®'* 

3.2260/R®*®®®® 

5.4491/R®*“®®® 

12.4377/R®***®®^ 






Table S. Boundary layer fchickneaa: Burger a flow 


e 

R * 2.0 

R - 2.828 

R « 3.464 

R ■ 4.0 

A(8)/R* 

ir 

0.3741 

0.4823 

0.4407 

0.4103 

O.8OO3/RO.‘»025 


0.3808. 

0.4878 

0.4434 DUA147 

0.8108/R®«‘»®'»1 

il”- 

.0.6018 

0.3043 

0.4603 

0.4283.. 

0.8427/R®***®®® 


0.6394 

0.5347 

0.4870 

0.4527 

0.9001/R®>^®2 


0.6988 

0.5823 

0.5291 

0.4912 

0.9911/R®‘5®®^ 

12’ 

0.7891 

0.6545 

0.5930 

0.5494 

1.1297/R®*5205 

!.• 

0.9269 

0.7646 

0.6902 

0.6381 

1.3422/R®*5368 

5 

12 

1.1444 

0.9385 

0.8437 

0.7780 

1.6Z86/r0.5552 

4 

12^ 

1.5101 

1. 2314 

ia023 

1.0138 

2.2435/R®.*^34 

3 

12^ 

2.1901 

1.7797 

1.5878 

1.4569 

3.2853/R®.5®67 

2 

12 

3.6943 

3.0067 

2.6809 

2.4587 

3.5407/R®*®®®® 

1 

TI" 

8.5786 

7.0497 

6.3159 

5.8136 

12.6387/R®*5®^® 
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is a substantial infirovsiiteat over the Oseen model in describing the flow 
near the cylinder. 


BURGBR8-£LO» FAST AN BbLXPTlC CYLINDER 


We consider first the case when- the-ma jot axis of the ellipse is 
paraUel.to the-uniform stream at infinity. The appropriate coordinate 
system is vhere 

x-^ cosh C cos n 


y ■ sinh K sin n 


(28) 


The curve 5 ■ Co represents an. ellipse with major axis 2 cosh Co along 
the x*>axis and. minor axis 2 sinh continuous inviscid potential 

flow which gives the convective velocity field satisfies the boundary value 
pro blem 

O 

T(C .n) - 0 

O 




y 


1 c 

T® 


sin n 


as C + * 


(29) 


The solution of equation (29) is 


Y(C*h) ■ e ® sinh (C -5 ) sin n 

o 


(30) 


If <i*(C,n) is the stream function for Burgers flow past an ellipse, 
the boundary value problem for tj* is defined as follows: 


I ^ n A'*' ^ 3’*' 3 \1 X 

—5. ♦ — ~ ♦ R I )}m(C,n) ■ 0 

ac^ \dc 3n 3n ac /) 


(31) 
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rip Tiff 


00 


(33) 


1 


i|)(5^,n) - -I (5-,n) • 0 

0 85 0 

i|» ~ — 0 ^ sin n 48 5 ♦ 

2 


The parameter R in equation (31) is proportional to the Reynolds nunber. 
Its value is 

R 

R 2 (34) 

cosh 5 ^ 

0 


is the Reynolds number based On the semimajor axis of the el- 
lipse. In the limit 5 ^ ^ 0, the ellipse degenerates to a finite flat 
plate. Since the potential flow parallel to a flat plate is a uniform 
stream* Burgers flow and Osaen flow are equivalent when 5 ^ ■ 0. 

As- in-the-case~jOf-the circle, the vorticity function has the form 


M (C »n ) 



R<KC,n) 



F (O T (n) 

n n n 


(35) 


where ^(C» 4 ) * e ° cosh (C * C ) cos n is the harmonic conjugate of ¥(C,n). 

o 

The resulting ordinary differential equations for F (O and T (h) are 

tt tt 

different versions of Mathieu's equation with solutions given by 


Pj^(0 • Cek, 




T (n) 


■ se, /n, - i- r 2 e 


(3b) 



» iiHiMMmt. W i n* ' 


If we make the following aubstltutlona In equations (39) and (36) i 

§ 

0«h» C«ike® (37) 

® 2 

the Burgers vorticity past an ellipso has exactly tho same form- as that post 
a circle: namelyi 

(.<..9 ) . -e® 8.k„ L - i .e„^9. - i (3») 

Poisson's equation [eq. (32)] becomes 

♦ LJL • - Jl h(s,6) o)(a,6) (39) 

3g2 302 2 

where 

h(z,6) « cosh 2(z 5^^) ** 20 (40) 

and its solution follows the technique described under "Solution." Omitting 
the details we find 

K It • 

4i(z ,0 ) “ e ® sinh e sin 0 4 JL / / h(Z|0) 

° ® ° 2 -s o 

• <u(z,0) G(e,0}z^,0^) dz d0 (41) 


where 


G(Zt0 *®0|0q) 


The coefficients in equation (3B) satisfy the infinite linear system 


(42) 
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V B r, ■ fi, , 

n«l 


k * li2)3|«tt 


where 


r 


kn 


^0 If *® 

e / / h<*,0)e"““ ^ ° C08h-B.c0fl 0 

ir 0 0 


• Oek„ 
n 




sin k0 se_ 
n 



- i c2 ) d*.., d0 
4 / 


(43) 


If the ellipse hes-its mSjor Sxis perpendicular to the flow at 
infinity, a different coordinate system is required. Consider the elliptic 
coordinates (m,X) defined by 


X ■■ sinh h cos.X 


y “ cosh y sin X (44) 

The curve y ** y^ represents an ellipse with major axis along the y~axis. 
The potential flow past this ellipse has stream function x(y>X) which 
satisfies exactly those conditions given in equation (2$). Ilius from 
equation (30) we have 

y 

x(y»X) ■ e ° sinh (y - y^) sin X (45) 


Since the potential flows in equations (30) and (43) are functionally the 
same, it follows that the Burgers solution in (y,X) coordinates will be 
similar to that in (C,n) coordinates. In fact, if the variables (y,X,x) 
are replaced by (C»n,¥), the problem is identical with that in equations 
(31) to (33), except that equation (32) is replaced by 


72^ 


1 0^t|) ^ \ 

sinh^C ♦ cos^n I 


w(f ,n) 


(46) 
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It follows that the atreaai futiotion for Burgers flow peat an ellipse whose 
major axis la perpendicular to the flow la given by equatlona (41) to (43) 
where h(8|d) Is new given by( 

h(»,0) ■ eofth 2(s ♦ 5^) ♦ 

The research that remains to bo done includes a dotaLLed analysis of 
the elliptical solutions. It is intention to coop^ito drag coefficients 
for various elliptical cylinders* examine separation phenomena, and calcu- 
late surface p.ressure distributions. The asymptotic form of the stream 
function will be used to deduce disp^lacement body shape and boundary Uyer 
thickness for various cases, the special cases of the flow- past a finite 
flat plate (both parallel to and perpendicular to the un'.form stream) will 
he Considered In detail and conpared with experime.ital data where possible. 
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Figure 1. 

StrecmllAOQ at the rear stagnation 
point P of the cylinder prior to 
separation 


V 





Figure 3. 

Convective velocity field at the 
rear stagnation point of the 
cylinder in Oseen flow 




Figure U. 

Convective velocity field in Burgers 
flow 





